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Motivated by recent experimental progress in the context of ultra-cold multi-color fermionic atoms 
in optical lattices, we have investigated the properties of the SU(A'’) Heisenberg chain with totally 
antisymmetric irreducible representations, the effective model of Mott phases with m < N particles 
per site. These models have been studied for arbitrary N and m with non-abelian bosonization 
[I. Affleck, Nuclear Physics B 265, 409 (1986); 305, 582 (1988)], leading to predictions about the 
nature of the ground state (gapped or critical) in most but not all cases. Using exact diagonalization 
and variational Monte-Carlo based on Gutzwiller projected fermionic wave functions, we have been 
able to verify these predictions for a representative number of cases with < 10 and m < N/2, and 
we have shown that the opening of a gap is associated to a spontaneous dimerization or trimerization 
depending on the value of m and N. We have also investigated the marginal cases where abelian 
bosonization did not lead to any prediction. In these cases, variational Monte-Carlo predicts that 
the ground state is critical with exponents consistent with conformal field theory. 


I. INTRODUCTION 

The possibility to load ultra cold fermionic atoms in 
optical lattices opens new perspectives in the investiga¬ 
tion of lattice models of strongly correlated systemsb^. 
When the optical lattice is deep enough, and when the 
fermionic atoms have an internal degree of freedom that 
can take N values (coming for instance from the nuclear 
spin in alkaline earths), the appropriate model takes the 
form of a generalized Hubbard model 


and to the fact that hifij is a constant in the manifold of 
states with one particle per site. 

More generally, when the number of degrees of free¬ 
dom is equal to iV >2, Mott phases can be realized for 
all integer values oim < N. The effective model now has 
SU(V) symmetry. This can be made explicit by intro¬ 
ducing the generators 
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which satisfy the SU(V) algebra 
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where dj^, and are creation and annihilation fermionic 
operators, hia = dlcOiai t is the hopping integral between 
pairs of nearest neighbors (f, j), and U is the on-site re¬ 
pulsion. 

When the average number of atoms per site m is an 
integer, and for large enough Ujt, the system is expected 
to be in a Mott insulating phaseb Fluctuations induced 
by the hopping term in the manifold of states with m 
fermions per site start at second order in tjU, and the 
processes that appear at this order consist in exchanging 
particles between pairs of neighboring sites, leading to 
the effective Hamiltonian 
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In the case of electrons with spin f or this Hamiltonian 
has SU(2) symmetry, and it is equivalent to the Heisen¬ 
berg model with coupling constant /U thanks to the 
identity 

E ®Lai/3a]/3«ia = 2^* • Sj + 

a,/3=t,4 


thanks to the identity 
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In the SU(V) language, working with m fermions per site 
corresponds to working with the totally antisymmetric 
irreducible representation (irrep) that can be represented 
by a Young tableau with m boxes in one column. For any 
allowed m, there is a conjugate equivalent representation: 
a system with m = k particles per site is equivalent to a 
system with m = N — k particles per site. 

The model Eq.2 captures the low-energy physics of 
multi-color ultra-cold atoms in optical lattices, systems 
for which remarkable progress has been recently achieved 
on the experimental side. For instance, the SU(A^)- 
symmetry has been observed in ultracold quantum gas 
of fermionic or ®^Sr''’. Another example*^ is the 

realization of one dimensional quantum wires of repul¬ 
sive fermions with a tunable number of components. 

The SU(A^) Heisenberg model with the fundamental 
representation at each site (m = 1), which corresponds 
to the Mott phase with one particle per site, has been 
investigated in considerable detail over the years. In one 
dimension, there is a Bethe ansatz solution for all values 
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oi N', and Quantum Monte Carlo simulations free of mi¬ 
nus sign problem have given access to the temperature 
dependence of correlation functions*’^’. In two dimen¬ 
sions, a number of simple lattices have been investigated 
for a few values of N with a combination of semiclassical, 
variational and numerical approaches, leading to a num¬ 
ber of interesting predictions at zero temperature 

In comparison, the case of higher antisymmetric irreps 
(to > 1) has been little investigated. In 2D, there is 
a mean-field prediction that chiral phases might appear 
for large to provided N/m is large enough^®’^^, and some 
cases of self-conjugate irreps such as the 6-dimensional 
irrep of SU(4) have been investigated with Quantum 
Monte Carlo simulations^®”'^' and variational Monte 
Carlo'®^. In ID, apart from a few specific cases'®^”"**^ , in¬ 
cluding more general irreps than simply the totally anti¬ 
symmetric ones, the most general results have been ob¬ 
tained by Affleck quite some time ago'^’**^. Applying 
non-abelian bosonization to the weak coupling limit of 
the SU(A^) Hubbard model, he identified two types of op¬ 
erators that could open a gap: Umklapp terms if > 2 
and N/m = 2, and higher-order operators with scaling 
dimension y = N{m — l)m~^ allowed by the Z^v/m sym¬ 
metry if N/m is an integer strictly larger than 2. This 
allowed him to make predictions in four cases: i) N/m 
is not an integer: the system should be gapless because 
there is no relevant operator that could open a gap; ii) 
N > 2 and N/m = 2: the system should be gapped 
because Umklapp terms are always relevant; hi) N/m is 
an integer strictly larger than 2 and y; < 2: the system 
should be gapped because there is a relevant operator al¬ 
lowed by symmetry. This case is only realized for SU(6) 
with TO = 2; iv) N/m is an integer strictly larger than 2 
and X > 2: the system should be gapless because there is 
no relevant operator allowed by symmetry. The only case 
where the renormalization group argument based on the 
scaling dimension of the operator does not lead to any 
prediction is the marginal case x = 2, which is realized 
for two pairs of parameters: (SU(8) to = 2) and (SU(9) 
TO = 3). These predictions are summarized in TAB. I. 
Finally, in all gapless cases, the system is expected to 
be in the SU(iV)j,^^ Wess-Zumino-Witten universality 
class with algebraic correlations that decay at long 
distance with a critical exponent rj = 2 — 2/N. 


1 Gapless 
? 

I Gapped 

TABLE I. Summary of the predictions of Refs. 41 and 42 
for a representative range of SU(A) with m particles per site. 
Note that models with m = k and m = N — k are equivalent 
up to a constant. Therefore the light gray shaded region can 
be deduced from the other cases and does not need to be 
studied. 

To make progress on the general problem of the SU {N) 
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Heisenberg model with higher antisymmetric irreps, it 
would be very useful to have flexible yet reliable numer¬ 
ical methods that would allow to test these predictions 
in a systematic way. In particular, the methods should 
not be limited to ID, or to cases such as self-conjugate ir¬ 
reps, for which there is a minus-sign free Quantum Monte 
Carlo algorithm. In this paper, we have decided to test 
the potential of Gutzwiller projected wave functions by 
a systematic investigation of the ID case discussed by 
Affleck using variational Monte Carlo (VMC). There are 
two main reasons for this choice. First of all, these wave 
functions have been shown to be remarkably accurate in 
the case of the SU(4) Heisenberg chain with the funda¬ 
mental representation by Wang and Viswanath'^'’’, who 
have in particular shown that they lead to the exact crit¬ 
ical exponent in that case. Besides, this approach can 
be easily generalized to higher dimensions, as already 
shown for the fundamental representation in a number 
of cases^'*’^"'. Moreover, it has been used by Paramekanti 
and Marston'®^ for the self-conjugate representation in 
one and two dimensions. 

In parallel, exact diagonalizations based on the exten¬ 
sion of a recent formulation by two of the present au- 
thors^^ will be used whenever possible to benchmark the 
VMC approach on small clusters and, in some cases, to 
actually conhrm the physics on the basis of a finite-size 
analysis. 

As we shall see, the combination of these approaches 
leads to results that agree with Affleck’s predictions 
whenever available, and to the identihcation of the sym¬ 
metry breaking pattern in the gapped phases. In addi¬ 
tion, it predicts that the marginal cases are gapless with 
algebraic correlations. 

The paper is organized as follows. The next section 
describes the methods, with emphasis on the variational 
wave functions that will be used throughout. The third 
section is devoted to a comparison of the results obtained 
using the simplest wave functions (with no symmetry 
breaking) with those of the Bethe ansatz solution for the 
TO = 1 case, with the conclusion that the agreement is 
truly remarkable. The fourth section deals with the cases 
where Umklapp processes are present {N > 2, N/m = 2), 
while the fifth one deals with the case where there is no 
Umklapp process but a relevant operator (SU(6) to = 2). 
The marginal cases are dealt with in the sixth section, 
and the case where N/m is an integer without relevant 
nor Umklapp operators in the seventh section. Finally, 
the critical exponents are computed and compared to 
theoretical values for all gapless systems in the eighth 
section. 


II. THE METHODS 

A. Gutzwiller projected wave functions 

The variational wave functions investigated in the 
present paper are obtained from fermionic wave functions 
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that have on average m particles per site by applying a 
Gutzwiller projector Pq that removes all states with a 
number of particles per site different from m: 

n 

i—1 p^m 

where n is the number of sites, and where the product 
over p runs over all values from 0 to iV except p = m. 

In the present paper, we will concentrate on simple 
fermionic wave-functions that, before projection, corre¬ 
spond to the ground state of trial Hamiltonians that 
contain only hopping terms. For SU(2), the inclusion 
of pairing terms have been shown to lead to significant 
improvements^®, but the generalization to SU(iV) is not 
obvious because one cannot make an SU (N) singlet with 
two sites as soon as iV > 2. In addition, in the case 
of the fundamental representation where Bethe ansatz 
results are available for comparison, these simple wave 
functions turn out to lead to extremely precise results as 
soon as IV > 2. 

In practice, the construction of a Gutzwiller projected 
wave function starts with the creation of a trial Hamilto¬ 
nian T that acts on n sites and is written with fermionic 
operators and aj^. When different colors are involved 
in T, and as long as there is no term mixing different col¬ 
ors, the Hamiltonian can be rewritten as a direct sum: 
T = Then, for each color, there will be one 

corresponding unitary matrix [/“ that diagonalizes Tq. 
So the new fermionic operators are given by: 

n n 

c- =y^C/“^d ='S^U°‘a^ 

and the trial Hamiltonian can be written in a diagonal 
basis: 

N n 

Oi—1 i—1 

with OJia < Wi+ia- 

In the Mott insulating phase, the system possesses 
nm/N particles of each color and exactly m particles per 
site. By filling the system with the nm/N lowest energy 
states of each color, the resulting fermionic wave function 
contains nm particles: 

N nm/N N mn/N n 

I'j^)=0 n "I io)=0 n E io) ( 4 ) 

Oi — 1 i—1 Q: = l i—1 j — 1 

in terms of which the variational wave function is given 
by 

|vI/g)=PS‘|^). (5) 

Since the Heisenberg model exchanges particles on 
neighboring sites, the simplest trial Hamiltonian that 


allows the hopping of particles and its corresponding 
Gutzwiller projected wave function are: 

n 

^Fermi ^ ^ + H.C.) ^ . 

i=l 

In cases where a relevant or Umklapp operator is present, 
the ground state is expected to be a singlet separated 
from the first excited state by a gap, and to undergo 
a symmetry breaking that leads to a unit cell that can 
accommodate a singlet. In practice, this means unit cells 
with d = N/m sites. To test for possible instabilities, we 
have thus used wave functions that are ground states of 
Hamiltonians that creates d-merization: 

n 

= X] + H.C.^ —)■ 1^^ (5)) . 

i^l 

Assuming that the mirror symmetry is preserved, the 
wave functions |'I'q(i5)) for dimerization {d = 2) and 
trimerization (d = 3) have only one allowed free parame¬ 
ter (5, and the hopping amplitudes in a unit cell are given 
by: 

ti = 1 — 5 ii i = d 

ti = l otherwise. 

To test for a possible tetramerization for SU(8) m = 
2, since the unit cell contains four sites, one additional 
free parameter is allowed (still assuming that the mirror 
symmetry is preserved in the ground state). Therefore, 
we have used the wavefunction | (di, ^ 2 )) with hopping 
amplitudes defined by: 

ti = 1 — Si if i = 2 
ti = 1 — S 2 if i = 4 
ti = 1 otherwise. 

This method is always well defined for periodic boundary 
conditions when N/m is even. But when N/m is odd, the 
ground state is degenerate for periodic boundary condi¬ 
tions if the translation symmetry is not explicitly broken, 
and one has to use anti-periodic boundary conditions for 
|v]/^ (Q)) (d = 2, 3) and |4/4 (0,0)). 

The hope is that if T is wisely chosen, then I'I'g) cap¬ 
tures correctly the physics of the ground state, i.e. with 
a good variational wave function, Eq = (Td H I'I'g) « 
Eq, the exact ground state energy. To check the perti¬ 
nence of this statement, we have compared the energies 
and nearest-neighbor correlations with those computed 
with ED on small systems with open boundary condi¬ 
tions. In the table TAB. H, one can see, for some sys¬ 
tems, the comparison between ED and VMG results for 
the ground state energy. The nearest-neighbor correla¬ 
tions will be compared in the next sections. Gonsidering 
the excellent agreement between the two methods for the 
cluster sizes available to ED, there are good reasons to 
hope that these Gutzwiller projected wave functions can 
quantitatively describe the properties of the ground state. 
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A 

m 

n 

ED 

VMC 

error [%] 

4 

2 

16 

-1.6971 

-1.6916 

-0.33 

4 

2 

18 

-1.6925 

-1.6866 

-0.35 

6 

2 

15 

-2.7351 

-2.7287 

-0.23 

6 

3 

12 

-4.0295 

-4.0261 

-0.08 

6 

3 

14 

-4.0162 

-4.0123 

-0.10 

8 

2 

12 

-3.1609 

-3.1587 

-0.07 

8 

2 

16 

-3.1857 

-3.1828 

-0.09 

9 

3 

9 

-6.0960 

-6.0810 

-0.25 

9 

3 

12 

-6.1162 

-6.0980 

-0.30 

10 

2 

15 

-3.3992 

-3.3919 

-0.21 


TABLE II. Comparison between the ED and VMC energies 
per site. The incertitudes on the VMC data are smaller than 
10~'*. The relative error is always smaller than 0.35%. 


B. Exact diagonalizations 

On a given cluster, the total Hilbert space grows very 
fast with N, and the standard approach that only takes 
advantage of the conservation of the color number is lim¬ 
ited to very small clusters for large N. Quite recently, two 
of the present authors have developed a simple method 
to work directly in a given irrep for the SU(7V) Heisen¬ 
berg model with the fundamental representation at each 
site^®, allowing to reach cluster sizes typical of SU(2) 
for any N. This method can be extended to the case 
of more complicated irreps at each site, in particular to¬ 
tally antisymmetric irreps, and the exact diagonalization 
results reported in this manuscript have been obtained 
along these lines. 


C. Correlation function and structure factor 


To characterize the ground state, it will prove useful 
to study the diagonal correlation defined by: 

a a 

( 6 ) 

The structure factor is then given by the Fourier trans¬ 
form of this function: 

C{k) = ^ --- y C (r) e*'^" (7) 

^ ^ 27r TO (iV - to) y ^ ^ 


where the prefactor has been chosen such that: 




n 

1'K 


(QMC) results. Since these systems are known to be 
gapless, I'E'q®’'™*) is the only relevant wave function to 
study. 

Let us start with the ground state energy. Using Bethe 
ansatz, Sutherland^ derived an exact formula for the 
ground state energy per site eo(iV) of the Hamiltonian 
Eq.2 that can be written as a series in powers of 1/A^: 


eo(iV) 


OO 

-1 + 2E 




(-i)^C(fc) 

Nk 


( 8 ) 


where ({k) = i® Riemann’s zeta function. 

eo{N) is depicted in FIG. 1 as a continuous line. The 
dashed lines are approximations obtained by truncating 
the exact solution at order N~^, k > 2. For compari¬ 
son, the variational energies obtained in the thermody¬ 
namic limit after extrapolation from finite size systems 
are shown as dots in FIG. 1. The agreement with the 
exact solution is excellent for all values of N, and it im¬ 
proves when N increases (see table TAB. HI). Quite re¬ 
markably, the variational estimate is better than the 
estimate even for SU(3). 


m = 1 



FIG. 1. Variational energy per site of SU(A) chains with 
the fundamental irrep at each site (dots) compared to Bethe 
ansatz exact results (solid line) and polynomial approxima¬ 
tions in 1/A (dashed lines). 


A 

BA 

VMC 

error [%] 

3 

-0.7032 

-0.7007 

-0.36 

4 

-0.8251 

-0.8234 

-0.21 

5 

-0.8847 

-0.8833 

-0.16 

6 

-0.9183 

-0.9173 

-0.11 

7 

-0.9391 

-0.9383 

-0.09 

8 

-0.9528 

-0.9522 

-0.06 

9 

-0.9624 

-0.9620 

-0.05 


III. SU(A) WITH m = 1 

In this section, we extend the SU (4) results of Wang 
and Vishwanath*® to arbitrary N for to = 1 (fundamen¬ 
tal representation), and we perform a systematic com¬ 
parison with Bethe ansatz and Quantum Monte-Carlo 


TABLE III. Comparison of the variational energies for m = 1 
systems obtained for infinite chains with exact Bethe ansatz. 
The incertitudes on the VMC data are smaller than 10“'*. 

We now turn to the diagonal correlations and its as¬ 
sociated structure factor defined by Eq.6 and Eq.7. At 
very low temperature, QMC has been used by Frischmuth 
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m = 1 n = 60 



FIG. 2. Comparison of the structure factors calculated 
with VMC (empty squares) and QMC (filled circles) for var¬ 
ious SU(A'^) systems. In the VMC calculations, anti-periodic 
boundary conditions have been used for SU(3) and SU(5), 
and periodic ones for SU(4). 


et aP for SU (4) and by Messio and Mila^ for various val¬ 
ues of N to compute this structure factor. The QMC 
data of Messio and Mila and the results obtained with 
VMC for n = 60 sites are shown in FIG. 2. Qualita¬ 
tively, the agreement is perfect: VMC reproduces the 
singularities typical of algebraically decaying long-range 
correlations. But even quantitatively the agreement is 
truly remarkable, and, as for the ground state energy, 
it improves when N increases. Clearly, Gutzwiller pro¬ 
jected wave functions capture the physics of the m = 1 
case very well. 


IV. SU(V) WITH m = N/2 

For these systems, there is a self-conjugate antisym¬ 
metric representation of SU(V) at each site. The ground 
states of such systems, referred to as extended valence 
bound solids'*^, are predicted to break the translational 
symmetry, to be two-fold degenerate and to exhibit 
dimerization since only two sites are needed to create 
a singlet, and the spectrum is expected to be gapped. 

We have investigated two representative cases, (SU(4) 
m = 2) and (SU(6) to = 3), with ED up to 18 and 14 
sites respectively, and the cases iV = 4 to 10 with VMC. 
The main results are summarized in FIG. 3. 

Let us start by discussing the ED results. Clusters 
with open boundary conditions have been used because 
they are technically simpler to handle with the method of 
Ref. 45, and because, in the case of spontaneous dimer¬ 
ization, they give directly access to one of the broken 
symmetry ground states if the number of sites is even. 
The gap as a function of the inverse size is plotted in the 
upper left panel of FIG. 3 for SU(4) and SU(6). In both 
cases, the results scale very smoothly, and a linear fit is 
consistent with a finite and large value of the gap in the 
thermodynamic limit. In the lower panels of FIG. 3, the 
bond energy is plotted as a function of the bond position 


N/m = 2 
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FIG. 3. ED and VMC results for various SU(iV) models 
with m = N/2. Upper left panel: size dependence of the en¬ 
ergy gap for SU(4) and SU(6). Upper right panel: optimal 
variational parameter 5 for periodic boundary conditions for 
SU(4), SU(6), SU(8), and SU(IO). Lower panels: energy per 
bond for SU(4) (left) and SU(6) (right) calculated with ED 
(circles) and VMC (squares) for open boundary conditions. 
Note that the optimal variational parameters 5opt are differ¬ 
ent in the upper right panel and in the lower panels because 
they correspond to different boundary conditions (periodic 
and open). 


for the largest available clusters (18 sites for SU(4), 14 
sites for SU(6)) with solid symbols. A very strong alter¬ 
nation between a strongly negative value and an almost 
vanishing (slightly positive) value with very little depen¬ 
dence on the bond position clearly demonstrates that the 
systems are indeed spontaneously dimerized. 

Let us now turn to the the VMG results. Since the 
relevant instability is a spontaneous dimerization, it is 
expected that the dimerized |'I'q( 5)) wave function al¬ 
lows one to reach lower energy than the one. 

This is indeed true for all cases we have investigated (up 
to V = 10 and to n > 100), and the optimal value of 
the dimerization parameter Jopt > 0 is nearly size inde¬ 
pendent and increases with N (see upper right panel of 
FIG. 3), in qualitative agreement with the gap increase 
between SU(4) and SU(6) observed in ED. To further 
benchmark the Gutzwiller projected wave functions for 
these cases, we have calculated the bond energy using the 
optimal value of 5 (open symbols in the lower panel of 
FIG. 3) for the same clusters as those used for ED with 
open boundary conditions. The results are in very good 
quantitative agreement. 

With the large sizes accessible with VMC, it is also 
interesting to calculate the diagonal structure factor de¬ 
fined in Eq.7. All the structure factors peak at fc = tt, 
but, unlike in the case of the fundamental representa¬ 
tion, there is no singularity but a smooth maximum (see 
FIG. 4). This shows that the antiferromagnetic corre- 
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lations revealed by the peak at fc = tt are only short 
ranged, and that the correlations decay exponentially at 
long distance, in agreement with the presence of a gap, 
and with the spontaneous dimerization. 

To summarize, ED and VMC results clearly support 
Affleck’s predictions that the N/m = 2 systems are 
gapped and point to a very strong spontaneous dimeriza¬ 
tion in agreement with previous results by Paramekanti 
and Marston'^^. 


N/m = 2 n = 120 



FIG. 4. Structure factor of various SU(A) models with 
m = N/2 calculated with VMC with the optimal variational 
parameter dopt. 


V. SU(6) WITH m = 2 

This case is a priori more challenging to study because 
the relevant operator that is generated in the renormal¬ 
ization group theory appears at higher order than the 
one-loop approximation. Therefore, the gap can be ex¬ 
pected to be significantly smaller than in the previous 
case. This trend is definitely confirmed by ED performed 
on clusters with up to 15 sites: the gap decreases quite 
steeply with the system size (see upper left panel of 
FIG. 5). It scales smoothly however, and a linear ex¬ 
trapolation points to a gap of the order AE ~ 0.2, much 
smaller than in the SU(6) case with m = 3 (AE ~ 4), 
but finite. On the largest available cluster, the bond en¬ 
ergy has a significant dependence on the bond position, 
with an alternance of two very negative bonds with a less 
negative one. 

These trends are confirmed and amplified by VMC. 
Indeed, the trimerized wave function IfE'g [5)') leads to a 
better energy for all sizes, and the optimal value scales 
very smoothly to a small but finite value i5opt ~ 0.03. 
This value is about an order of magnitude smaller than 
in the SU(6) case with m = 3, but the fact that it does 
not change with the size beyond 60 sites is a very strong 
indication that the system trimerizes (by contrast to the 
marginal case shown in FIG. 9). The trimerization is con¬ 
firmed by the lower plots. For n = 15, the VMC results 
are again in nearly perfect agreement with ED, and for 
n = 60, the bond energy shows a very clear trimerization. 

To test the nature of the long-range correlations is of 


SU(6) m = 2 



i 


FIG. 5. ED and VMC results for the SU(6) model with 
m = 2. Upper left panel: size dependence of the energy 
gap. Upper right panel: optimal variational parameter 5 for 
periodic boundary conditions. Lower left panel: energy per 
bond calculated with ED (circles) and VMC (squares) on 15 
sites with open boundary conditions. Note that the optimal 
variational parameter 5 = 0 in that case. Lower right panel: 
energy per bond calculated with VMC with periodic boundary 
conditions. 


course more challenging than in the previous case since a 
small gap implies a long correlation length. And indeed, 
on small to intermediate sizes, the structure factor has a 
sharp peak at fc = 27r/3 very similar to the SU(3), m = I 
case. However, going to very large system sizes (up to 



0 27r/3 TT 1.957r/3 2.257r/3 

k k 



FIG. 6. Upper left panel: Structure factor of the SU(6) 
model with m = 2 calculated with VMC using a trimerized 
wave function with the optimal variational parameter. Upper 
right panel: zoom on the region near k = 27r/3. It clearly 
shows that the structure factor is smooth. Lower panel: real- 
space diagonal correlations for 60 sites. 
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n = 450 sites), it is clear that the concavity changes sign 
upon approaching k = 27r/3 (see upper right panel of 
FIG. 6), consistent with a smooth peak, hence with ex¬ 
ponentially decaying correlation functions (see also lower 
panel of FIG. 6). 

In that case, in view of the small magnitude of the gap, 
hence of the very large value of the correlation length, 
it would be difficult to conclude that the system is defi¬ 
nitely trimerized on the basis of ED only. In that respect, 
the VMC results are very useful. On small clusters, the 
Gutzwiller projected wave function with trimerization is 
nearly exact, and VMC simulations on very large systems 
strongly support the presence of a trimerization and of 
exponentially decaying correlations"^®. 


VI. MARGINAL CASES: SU(8) WITH m = 2 AND 
SU(9) WITH m = 3 

These two systems are the only ones which possess op¬ 
erators with scaling dimension y = 2. They are therefore 
the only cases where it is impossible to predict whether 
the system is algebraic or gapped on the basis of Affleck’s 
analysis. As far as numerics is concerned, these cases can 
again be expected to require large system sizes to con¬ 
clude. 






FIG. 7. ED and VMC results for the marginal cases SU(8) 
with m = 2 and SU(9) with m = 3. Upper left panel: size de¬ 
pendence of the energy gap for both cases. Upper right panel: 
optimal variational parameter 5 for the SU(9) case with pe¬ 
riodic boundary conditions. The results for SU(8) are not 
shown because they identically vanish for periodic boundary 
conditions. Lower left panel: energy per bond for SU(8) cal¬ 
culated with ED (circles) and VMC (squares) for open bound¬ 
ary conditions. Note that the optimal variational parameters 
5opt are different from zero with open boundary conditions. 
Lower right panel: energy per bond for SU (9) calculated with 
ED (circles) and VMC (squares) for open boundary condi¬ 
tions. 


The ED results are quite similar to the previous case. 
The scaling of the gap is less conclusive because the last 


three points build a curve that is still concave and not 
linear like in the previous case (see the upper right panel 
of FIG. 7). So one can only conclude that if there is a 
gap, it is very small, especially for SU(8) with m = 2. 
The bond energies build a pattern which is consistent 
with a weak tetramerization for SU(8) with m = 2, and 
with a significant trimerization comparable to the SU (6), 
m = 2 case for SU(9) m = 3. 

The VMG method turns out to give a rather different 
picture however. For SU(8) with m = 2, two variational 
wave functions (Ji, ^ 2 ))) can be tested. 

Interestingly, for n = 16 with open boundary conditions, 
I^Fermi^ fails to reproduce the bound energies pattern 
observed with ED but (0.054, —0.036)) is successful 
(see lower left pannel of FIG. 7). This pattern, which 
could be interpreted as a weak tetramerization, is in fact 
probably just a consequence of the four-fold periodicity 
of algebraic correlations in the presence of open bound¬ 
ary conditions. Indeed, it turns out that, for any system 
size with periodic boundary conditions, the minimization 
of the energy using (^i,(52)) failed to hnd a solution 
for any |(5i_2| > 0.002. Therefore is believed to 

be the best variational wavefunction. The conclusion is 
that there is no tetramerization, and that the correlations 
must be algebraic. This is also supported by the struc¬ 
ture factor, which seems to have a singularity at A: = 7r/2 
(see FIG. 8). 
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FIG. 8. Upper panels: Structure factor of the SU(8) model 
with m = 2 (left) and of the SU(9) model with m = 3 
(right) calculated with VMC with periodic boundary condi¬ 
tions. Lower panels: real space correlations. The four plots 
represents results obtained with 


Let us now turn to SU(9) with m = 3. This system 
could in principle be trimerized, and therefore |'Fq’'“‘) 
and |4 'q(( 5)) have been compared. For small clusters, 
there is a large optimal value of <5, actually much larger 
than for SU(6) with m = 2, and the bond energies are 
typical of a strongly trimerized system, in agreement with 





















ED. However, Jopt decreases very fast with n until it van¬ 
ishes for n > 100 whereas, for SU( 6 ) with m = 2, Jopt 
levels off at a finite value beyond n = 60 (see FIG. 9). 
We interpret this behavior as indicating the presence of 
a cross-over: on small length scales, the system is effec¬ 
tively trimerized, but this is only a short-range effect, and 
the system is in fact gapless with, at long-length scale, 
algebraic correlations. 

SU(6) m = 2 SU(9) m, = 3 
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FIG. 9. Energy per site as a function of the variational 
parameter 5 for the SU( 6 ) with m — 2 (left) and SU(9) with 
m = 3 (right). 

One can again calculate the structure factor using the 
best variational wave function (in both cases 
for big enough systems) to check if a discontinuity exists. 
The results displayed in the upper plots of FIG. 8 clearly 
show a discontinuity at fc = 7 r /2 for the SU( 8 ) and at 
k = 27r/3 for SU(9). These discontinuities indicate an 
algebraic decay of the long-range correlations. The lower 
plot shows that even if these systems are gapless, there 
is a maxima of the correlation every N/m sites. 
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VII. EXAMPLE WITH IRRELEVANT 
OPERATOR: SU(IO) WITH m = 2 

For completeness, we have also looked at a case where 
there is an irrelevant operator of scaling dimension larger 
than 2, namely SU(IO) with m = 2. As expected, the 
best variational wave function is for all sizes, 

and the structure factor exhibits discontinuities at k = 
27r/5, consistent with a gapless spectrum and algebraic 
correlations. 


VHI. CRITICAL EXPONENTS 


predicts an algebraic decay of the long-range correlations 
function according to: 

n ( \ , Cfc cos (2Trrm/N) 

- ;:2 ^ 

where rj = 2 — 2/N is the critical exponent. 

For systems with periodic boundary conditions, one 
can define two distances between two points, which nat¬ 
urally leads to the following fitting function*^: 

co(r“‘" + {n — r)~'')+Ck cos {2'Krm/N) (r”'* + {n — r)~^) 

with four free parameters: cq , v and Ck, rj, the amplitudes 
and critical exponents of the components at A: = 0 and 
k = 2'Km/N respectively. 

There is a large degree of freedom in the choice of the 
fitting range. One could in principle select any arbitrary 
range of sites [xi,Xf], 0 < Xi < Xf < n—1. The problem 
is that each range will give different critical exponents. 
In order to obtain some meaningful results the following 
method has been chosen. Using the periodicity of the 
systems, only the ranges with Xi = a and a;/ = n — a — 1, 
1 < a < n/2, have been considered. For each value of a, 
the coefficient of determination of the fit has been com¬ 
puted and if its value is higher than 0.999 than the range 
[a, n — a — 1] is selected to perform the extrapolation of 
the critical exponents. If the value is too low, the fit is 
considered to be bad and the range with a •(— a -I- 1 is 
tested. If no good range can be found with this crite¬ 
rion, the condition over the coefficient of determination 
is relaxed to be higher than 0.995 and the first fit with a 
residual sum of squares divided by n that is smaller than 
10 “^ is selected. 
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FIG. 10. Gritical exponents rj of the gapless systems as a 
function of the system size. The squares, circles and triangles 
correspond respectively to m = 1,2, and 3 particles per site. 
All values given here have been calculated with 


Motivated by the remarkably accurate results ob¬ 
tained in previous works for the case the fundamental 
representationwe have tried to use the VMG results 
to determine the critical exponent that controls the decay 
of the correlation function at long distance, Eq.7. For the 
particular case of gapless systems, conformal field theory 


The critical exponents rj obtained in this way are shown 
in FIG. 10. The theoretical values of the critical expo¬ 
nents 77 = 2 — 2/N are shown as straight lines. In all 
cases, the extracted exponents agree quite well with the 
theoretical predictions when n is large enough. In partic¬ 
ular, for a given V, the exponent rj does not depend on 
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TO, as predicted by non-abelian bosonization. The criti¬ 
cal exponents v has also been extracted but, as already 
observed®, a precise estimate is difficult to get. Never¬ 
theless, for TV = 3,4, z/ S [1.8,2.25] and for N > 5, 
V G [1.95,2.05] for the largest systems. 

IX. CONCLUSIONS 

Using variational Monte Carlo based on Gutzwiller 
projected wave functions, we have explored the prop¬ 
erties of SU(iV) Heisenberg chains with various totally 
antisymmetric irreps at each site. In the case of the fun¬ 
damental representation, which is completely understood 
thanks to Bethe ansatz and to QMC simulations, these 
wave functions are remarkably accurate both regarding 
the energy and the long-range correlations. In the case 
of higher antisymmetric irreps, where field theory argu¬ 
ments are in most cases able to predict that the system 
should be gapless or gapped, allowing for a symmetry 
breaking term in the tight binding Hamiltonian used to 
define the unprojected wave function leads to results in 
perfect agreement with these predictions, and the ground 
state is found to be spontaneously dimerized or trimer- 


ized. Finally, in the two cases where the operator that 
could open a gap is marginal, SU(8) with to = 2 and 
SU(9) with TO = 3, this variational approach predicts 
that there is no spontaneous symmetry breaking, and 
that correlations decay algebraically. These results sug¬ 
gest that the operators are marginally irrelevant in both 
cases. It would be interesting to test these predictions 
either analytically by pushing the renormalization group 
calculations to higher order, or numerically with alterna¬ 
tive approaches such as DMRG or QMC. 

In any case, these results prove that Gutzwiller pro¬ 
jected fermionic wave functions do a remarkably good 
job at capturing quantum fluctuations in one-dimensional 
S\]{N) Heisenberg models with totally antisymmetric ir¬ 
reps. Gonsidering the encouraging results obtained in 2D 
for the SU(A^) Heisenberg model with the fundamental 
irrep at each site on several lattices, one can legitimately 
hope these wave functions to be also good for the SU(TV) 
Heisenberg model with totally antisymmetric irreps at 
each site in 2D. Work is in progress along these lines. 
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